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ABSTRACT 

We calculate the full stellar-structural evolution of donors in AM CVn systems formed 
through the WD channel coupled to the binary's evolution. Contrary to assumptions 
made in prior modelling, these donors are not fully convective over much of the AM 
CVn phase and do not evolve adiabatically under mass loss indefinitely. Instead, we 
identify three distinct phases of evolution: a mass transfer turn-on phase (during which 
-forb continues to decrease even after contact, the donor contracts, and the mass trans- 
fer rate accelerates to its maximum) , a phase in which the donor expands adiabatically 
in response to mass loss, and a cooling phase beginning at P or b ~ 45-55 minutes dur- 
ing which the donor contracts. The physics that determines the behaviour in the first 
and third phases, both of which are new outcomes of this study, are discussed in some 
detail. We find the overall duration of the turn-on phase to be between ~ 10 4 -10 6 yrs, 
significantly longer than prior estimates. We predict the donor's luminosity, L, and 
effective temperature, T e g. During the adiabatic expansion phase (ignoring irradiation 
effects), L w 1O" 6 -1O" 4 L and T cS « 1000-1800 K. However, the flux generated in 
the accretion flow dominates the donor's intrinsic light at all times. The impact of ir- 
radiation on the donor extends the phase of adiabatic expansion to longer P or b , slows 
the contraction during the cooling phase, and alters the donor's observational charac- 
teristics. Irradiated donors during the adiabatic phase can attain surface luminosities 
up to « 10~ 2 Lq. We argue that the turn-on and cooling phases both will leave signif- 
icant imprints on the AM CVn population's P or t,-distribution. Finally, we show that 
the eclipsing AM CVn system SDSS J0926+3624 provides evidence that WD-channel 
systems with non-zero entropy donors contribute to the AM CVn population, and we 
discuss the observational signature of the donor in this system. 

Key words: binaries: close — gravitational waves — stars:individual (RX J0806+1527, 
RX J1914+2456, SDSS J0926+3624)— white dwarfs 
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mass, essentially pure-He donor onto a white dwarf (WD) 
accretor is taking place. Thus, AM CVn systems form the 



ries initiating mass transfer. The characteristics of the AM 

CVn population can thus provide insights into the physics 
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for strong constraints from the AM CVn population on our 
theories of binary formation and evolution are, thus, very 
bright. In the meantime, advances from electromagnetic ob- 
servations and developing theory will begin addressing these 
same questions. 

There have been three distinct binary evolution chan- 
nels proposed to form AM CVn systems. Two of these forma- 
tion channels involve a series of two common-envelope (CE) 
events which bring the remnant cores of a main-sequence 
(MS) binary close enough for GW-emission to drive them 
into contact, initiating the AM CVn phase. These channels 
are disti nguished by the state of the proto-donor prior to 
contact ()Nelemans et al.ll200ll ). In what we'll call the He- 
star channel, the proto-donor is able to ignite He prior to 
contact. In the other, which we'll refer to as t he WD chan- 
nel, t he donor is partially-to-very degenerate (|Delove et al.l 
l2005l ) and has not undergone any He burning. The third 
channel involves an evolved MS star with a WD compan- 
ion starting stable m ass transfer just before it evolves up 
the red-giant branch (|Podsiadlowski et al.ll20olj| ). The star's 
core is dominated by He, allowing the binary to evolve to ul- 
tracompact configurations. We'll refer to this third channel 
as the evolved-MS channel. 

In mass transferring binaries, the donor's structural 
response to mass loss plays a central role in determin- 
ing t he binary's evolution (see, e.g., iDelove fc Bildstenl 
120031 ). The extent to which the donors from each forma- 
tion channel have been modelled varies considerably. He- 
star channel donors have only been modelled by a fit 
l|Nelemans et al.l l200ll) to a single relev a nt bi nary evolu- 
tion calculation by iTutukov fc Fedoroval (|l989l ). The WD 
channel donors have been considered more extensively. 
iNelemans et al.l (|200ll ) modelled these donors as fully- 
degenerate WDs obeying the zer o-temperature He WD 
mass-radius (M-R ) relation. Later. [Delove et al.l (|2005T ). us- 
ing donor structural models developed for ultr acompact low- 
mass X-ray binaries (jDelove fc Bild stcn 2003!) , modelled the 
evolution of non-zero entro py donors from the WD-channel. 
iPodsiadlowski et al.l (|2003h has presented a detailed study 
of the evolved-MS star channel. 

In this paper, we focus on the WD channel systems. 
The lDelove et all (|2005l ) study had several limitations: they 
assumed the donors to be fully convective and to evolve adia- 
batically under mass loss due to the very large mass transfer 
rates, M, produced in AM CVns. With these assumptions, 
the evolution is completely determined by the total mass of 
the two components and the donor's specific entropy, which 
determines the M-R relation the donor follows and, amongst 
other system parameters, the binary's M(P or b) evolution. 
The inclusion of non-zero entropy donors in this population 
resulted in shifting systems to larger M at fixed P or b and 
allowed systems to evolve out to longer P or b as compared 
to a population of zero-temperature donors. This provided 
an observational diagnostic of an AM CVn system's initial 
conditions. The IDelove et al.l (|2005l ) calculation also indi- 
cated that WD channel systems with hot donors and He- 
star channel systems would not be distinguishable based on 
M determinations alone. 

Since the entropy structure of the donor determines 
its response to mass loss (see the discussion in M4. 1 ■ 1 1 and 
Appendix [A")l . assuming a fully convective structure will 
overestimate the binary's P or b evolution rate if the donor's 



true structure in not adiabatic. Further, while assuming the 
donor responds adiabatically to mass-loss is certainly valid 
at early stages of the AM CVn phase when M is very large, 
this assumption may break down as the binary evolves out- 
ward in P or b to lower M. 

To rectify these shortcomings and provide predictions 
for the donor's luminosity and effective temperature, we un- 
dertook a study to model the donors within the context of 
a full stellar structural calculation coupled to the evolution 
of the AM CVn binary. In particular, we sought to deter- 
mine the donor's observational signatures as a function of 
initial system parameters and P or b, as well as to determine 
how a complete treatment of the donor's structure alters the 
conclusions of earlier studies. We also sought to model accu- 
rately the critical phase of mass transfer turn-on, when the 
donor is coming into contact and a transition in the P or b- 
evolution from GW-dominated in-spiral to the mass-transfer 
dominated expansion occurs. This phase of WD channel sys- 
tem evolution has not been modelled in any detail prior to 
this study. 

To carry out the computations, we developed a new 
stellar evolution code to perform the coupled donor/binary 
calculations. A brief description of this code, detailing our 
numerical modelling of the donor and binary is provided 
in fJ2] In 331 we describe the procedure for determining the 
range of initial donor degeneracy and donor structural mod- 
els used for our AM CVn phase calculations. We present the 
results of our calculations in i|3]and discuss these results and 
several of their applications in ij5] Finally, we summarize in 

m 



2 DETAILS OF OUR NUMERIC 
CALCULATIONS 

To achieve a necessary level of flexibility in terms of input 
physics and defining system of equations for this and related 
projects, we implemented a stellar evolution code capable of 
handling most ID stellar evolution calculations expressed as 
boundary-value problems. This code is written in C++ and 
is comprised of a suite of abstract classes defining the nec- 
essary components of a stellar evolution calculation. At its 
core is a generalized ID relaxation-meth od boundary-value 
problem (BVP) solver (see, e.g, §17.3 of lPress et _aLlll992D 
that utilizes the sparse-matrix solver UMFPACK0 to per- 
form the necessary matrix inversions. The defining system 
of ordinary differential equations (ODEs) and necessary in- 
put physics are implemented as separate abstract classes, 
providing the desired flexibility. Utilizing this new code, we 
solve for both the donor's structure and binary parameters 
within a single set of relaxation iterations at each time step. 
Below we detail our numerical treatment for both compo- 
nents of our system and describe the input physics. 

2.1 Calculation of the Donor's Structure 

For the donor, we solve the standard set of ID hydrostatic 
stellar structure equations augmented by an automatic mesh 

2 http:/ /www. cisc.ufi.edu/research/sparse/umfpack/ 
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allocation algorithm l|Eggletonl Il97lh on a 200-point non- 
lagrangian mesh: 
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Here p, T, P, r, and are the density, temperature, pressure, 
radius, and luminosity at each mesh point, m is the mass 
interior to each mesh point, t is time, e is the local nuclear 
energy generation rate (which we set to zero throughout 
this work), and cp is the specific heat at constant pressure. 
In equation <(Sj — which determines the mesh allocation — 
P c and P e are the pressures at the stellar centre and exte- 
rior point, respectively, ai, a.2 are constants that weight the 
importance of pressure and mass gradients in determining 
the mesh spacing, and 9 is an overall normalization con- 
stant determined implicitly at each time step. The quan- 
tity V a d = (d\nT/d\nP) s , where s is the specific entropy, 
and V = d In T/d InP describes the actual run of T with 
P within the donor. We take V = V ra d when V ra d < V a d, 
where 



V ra d = 



Pk I 



IG-KacG T 4 m 



(6) 



and k is the Rosseland mean opacity. Otherwise we deter- 
mine V using mixing length theory with the mixing length 
set equal to one pressure scale height. 

The central boundary conditions are determined by the 
standard central expansions of the stellar structure equa- 
tions. At the outer boundary, we equate the values of the it- 
erated quantities at the outermost mesh point, T e , P e , r e , l e , 
and m e with the corresponding values at the photosphere 
(i.e., we do not calculate a separate atmosphere model). 
Grey atmospheres and the Eddington approximation are as- 
sumed. The photospheric pressure is therefore given by: 



phot q 



2 GM 2 



3P? 

^pho 



(T) 



where K p hot = «(ft>hot, T e e), p P hot is the photospheric den- 
sity, and Teg is the donor's effective temperature. This gives 
the matching conditions T e — T e g, P e = P p h ot , r e = R2, and 
m e = M2, and l e = L. Here M2, R2 and L are the donor's 
total mass, radius, and luminosity, with L = A-KaR^T^. 

To produce a smoothly varying outer boundary condi- 
tion, equation J7J is not solved directly each time photo- 
spheric values are required. Instead, a 3 x 3 grid of (T e g,g — 
GM2/R2) points is constructed that straddle the donor's 
current photospheric conditions. We calculate P prio t from 
equation (0 at these values and use these to determine 
Pphot by cubic-spline interpolations during the iterations. If 
a point outside of the current grid is needed, the grid is re- 
centred and new P p hot values calculated. The T c g and g grids 
have spacings Alog(T eff /K) = 0.05, Alog(p/cms- 2 ) = 0.05. 



2.2 Calculation of the Binary Parameters 

To specify the binary system, we add the following three 
equations to the above set of ODEs: 



d In M 2 
dt 



M2 
M 2 



M 

'Ma 



d In Mi _ Mi 

Jt ~~ mT ' 



(8) 
(9) 



and 



d In a ^ J 

dt ~ J 



Mi 

mT 



Mi + Ah 
(Mi + 2M 2 



+ 7 /(M 2 +2Mi) 
JV12 



(10) 



where Mi is the accretors mass, a is the orbital separa- 
tion, J is the orbital angular momentum. The quantities 
Mi, M 2 specify the mass evolution rate of the accretor and 
donor; both are, in general, functions of the binary's and, 
possibly, the donor's state. We express the a evolution in 
the form of equation (|10|l to allow for this generality. How- 
ever, here, we consider conservative mass transfer, setting 
Mi = M = —M2 always. We con sider only gravity wave 
emission (|Landau fc Lifshitj Il97ll ) contributions to the J 
evolution, ignoring the possibility of accretor spin-up in the 
J evolution. This latter effect may play a significant role in 
determinin g whether an AM CVn binary survives th e onset 
of contact (jNelemans et al.l l200ll : iMarsh et al.ll2004T ). This 
omission is justified since we are interested in understand- 
ing the parameter space potentially allowed to AM CVn 
systems, not in a detailed characterization of the actual, re- 
alized population. This latter question will be the subject of 
future studies. 

T o dete rmine M2 we u t ilize t he prescriptions of iRitteJ 
1 19881 ) and iKolb fc Ritteil (|l990l ). Specifically, when the 
dono r's Roche radius, Rl, here calculated as (|Paczvnskil 
[1963) 



Rl « 0.46a 



M 2 



Mi + M 2 



1/3 



(11) 



(valid for q = M2/M1 < 0.8) is larger than R2, the mass 
loss is modelled as an isothermal flow of a classical gas with 
a rate 



M 2 = Mo exp 



Z AP 
Up 



(12) 



where Mo > is the mass transfer rate off the donor 
when R2 = Rl and depends on the surface properties 
(|Ritteilll988l ). Hp is the pressure scale height at the pho- 
tosphere, and AP = R2 — Rl- When Rl < R2, we model 
the sub-photospheric contributions with an adiabatic flow 
jKolb fc Ritterlll99(t). T o do so, we rewrite equation (A17) 
of lKolb fc Ritteil (|l990h in a more general form to allow for 
degeneracy or non-ideality: 



-M 2 = M +27rPi(l/g) 



Rl 
GM 2 



(B2-AB) 



(13) 



where Pi and F2 are defined in IKolb fc Ritteil (|l990l ) and 
Ti = (91nP/<91np) s . We determine M implicitly at each 
timcst ep, a method that great ly improves the numerical sta- 
bility (|Buning fc Ritteill2006l ). 
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2.3 Input Physics 

For the donor's composition, we assume a fixed, homoge- 
neous mixture with He, C, N, and O mass fractions of 
0.981, 1.41 x 1(T 4 , 0.0122, and 1.06 x 10" 3 (typical of the 
core composition of solar metallicity stars at the base of 
the red-giant branch that have undergone C NO-cycIe burn- 
ing; [SchSeLelalE^ The remaining 
metals are assumed to have a solar abundance pattern . We 
use OPAL radiative opacities (|lglesias fc Rogers! [l996h cal- 
culated for this mixture for temperatur es above 10 4 K. Be- 
low 10 4 K, we use iFerguson et~ail (|2005h opacities calculated 
excluding contributions due to grains (which form beginning 
at T ~ 2500 K ). Where conductive o pacities are relevant, we 
utilize those of lPotekhin et all |l999h . A s metals make negli- 
gible contributions to the donor's EOS (IChabrier fc Baraffg 
1 19971 ). the pure-He EOS of lWinisdoerffer fc Chabrieri (|2005l ) 
supplemented by radiation contributions are used. This EOS 
is calculated using a free-energy minimization model and 
provides the thermodynamic quantities of a He fluid over a 
wide range of density and temperature including, in partic- 
ular, the regimes where He undergoes thermal and pressure 
ionization state transitions. 



3 THE INITIAL STATE OF THE DONOR IN 
WHITE DWARF CHANNEL SYSTEMS 



As shown in lDelove et al.l (|2005h . evolution during the post- 
contact AM CVn phase is strongly influenced by the state 
of the donor at contact. In this section, we determine a 
set of initial donor models that encompass the range of 
donor properties at contact produced in WD channel AM 
CVn systems and discuss the expected distribution of these 
properties based on th e population synthesis calculation of 
INelemans et al.l (|200ll V 

3.1 Determination of the Initial Donor Models 

To determine the range of parameters our in itial donor mod- 
els sh ould cover, we utilize data from the INelemans et al.l 
ll200l[) population s ynthesis model and the methodology of 
iDelove et al.l (|2005l ). Specifically, for each WD channel sys- 
tem in the population synt hesis model that is exp ected to 
survive initial contact (see INelemans et al.ll200lf ), we de- 
termine the proto-donor's core conditions at the beginning 
of the second CE using existing ste llar evo l ution calcula- 
tions performed w ith the EZ code (|Paxtonl 120041 ') for the 
IDelove et al.l (|2005l ) study. During the rapid CE phase, the 
proto-donor's central degeneracy doesn't change. So we uti- 
lize the central degeneracy parameter — ip c = EF, c /kT c ~ 
p c /(1.2 x 10~ 8 T 3/2 )K 3/2 cm 3 g -1 , where p c , T c , E F , C are the 
central density, temperature, and electron Fermi energy — 
to map between the EZ stellar models pre-CE and a set of 
post-CE isolated He WD models calculated using with our 
own stellar evolution code. 

Once we determine this post-CE donor model, for each 
system we evolve the donor's radius and the orbital separa- 
tion from this post-CE state to determine when GW emis- 
sion drives the binary back into contact. The population 
synthesis data provides the post-CE orbital separation, ao, 
and the initial accretor and donor masses, Mn and Mi i. 



We use these quantities and our single He-star tracks to 
calculate, self-consistently, the time elapsed in the post-CE 
phase, fcontact, before R2 = Rl- This determines the donor's 
contact radius and i/j c (which we denote by J?2,i and V'c.O- 
When taken over the set of population synthesis data, this 
procedure allows us to determine the range of ip c ,i, M2,i, and 
Mi,i that occurs in the WD channel systems. 

Before discussing the distribution of donor contact pa- 
rameters obtained with this procedure, we mention one com- 
plication we encountered. In some systems, the standard 
CE-evolution prescription (see, e.g. JWeb bink 1984]) predicts 
post-CE conditions where the donor is already in Roche con- 
tact. This is due to some systems having sufficiently hot 
proto-donors with R2 significantly greater than the zero- 
temperature configuration. It is unclear whether such sys- 
tems will cleanly exit the CE or simply be driven to a prompt 
merger. 

A realistic determination of this question is beyond the 
scope of this paper. Instead, in order to examine how this is- 
sue could influence the distribution of donor contact param- 
eters, we considered two different, rough criteria for deter- 
mining which systems survive the CE. In one, we removed all 
such systems from the population. In the second, we looked 
at the evolution of R2 (based on the isolated He WD struc- 
tures) and Rl to determine if 7?2 < Rl was ever satisfied 
before the system was driven to a — due to GW emission. 
This could happen if the proto-donor is able to contract 
more rapidly than GW emission decreases Rl- 

We found that the choice of criteria effected mainly 
the distribution of tp c ,i at fixed Af2,i for systems with 
M 2 ,i < 0.175M©. Taking the second criteria, About 50% 
of systems with 0.10M© < M 2 ,< < 0.125M Q and 25% of 
those with 0.125M© < M 2>i < 0.175M© do not exit the 
CE by this rough criteria, where as only 4% of those with 
M2,i > O.175M0 are affected. More importantly, the range 
of tp c ,i in the overall population is not affected by these 
border line systems, since donors with M2,i > O.2Af0 ex- 
hibit the full range of -0 c ,i and these more massive donors 
all appear to exit the CE cleanly. This general conclusion 
also applies to concerns about donors with M2,i ~ O.IMq. 
The progenitors of these donors will not have developed 
a distinct core/envelope structure at the start of the CE, 
making it unlikely that they will survive the CE-event 
l|Taam fc SandquistJ |2000T ). Thus, excluding such systems 
from consideration will also not influence the 4> c ,i range ex- 
pected in this population. 

3.2 The Distribution of Initial Donor Parameters 

The results of the above exercise produces a population 
with M 2 ,i « 0.1 - 0.325M Q and log(^ c ,i) ~ 1.0 - 4.0 (i.e., 
donors are mildly to extremely degenerate). Our determina- 
tion of the donor's p c , T c at contact for each system in the 
INelemans et al.l (|200ll ) population synthesis model, is shown 
in Figure [T] by crosses. Each cross represents the starting 
point for a number of systems, so the density of crosses does 
not correspond to the number of systems starting within 
a given region of this parameter space. Dashed lines show 
lines of constant ip c , while solid curves show the evolution 
along our isolated He WD tracks. We have excluded sys- 
tems not exiting the CE according to our second criteria 
above from this plot; if we were to use the first criteria, the 
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log(p c ) 

Figure 1. The don ors' p c , T c at contact for each WD channel AM 
CVn system in the iNelemans et al.l ll200ll ) population synthesis 
model (crosses). For comparison, the solid lines show the p c , T c 
evolution along each of our isolated He-star evolution tracks (with 
M 2 values of, from bottom to top, 0.10, 0.15, 0.20, 0.25, 0.30, and 
0.325 M ). Lines of constant logi/> c =1.0, 2.0, 3.0, and 4.0 are 
shown by the dashed lines (top to bottom). The solid circles show 
the central conditions of the set of initial donor models used to 
begin our AM CVn phase calculations. 




5 10 15 20 

P orb (minutes) 

Figure 2. The contact P or b-distribution as determined by self- 
consistent treatment of the donor's cooling before contact (solid 
histogram) compared to assuming all donors are fully degenerate 
at contact (dash dotted histogram). The number of systems in 
each bin gives the total number of cu rrent WD channel AM CVn 
systems in the lNelemans et al.l J200lj) population synthesis model 
whose contact P or b fell within each bin. The solid and dashed 
curves show the corresponding cumulative fractional distribution 
of contact P or b- 



only change here would be the loss of most systems with 
0.1 < M 2 ,i < 0.15M Q and 1.0 < log(V>c,i) < 1-5. The solid 
circles in this figure show the central conditions of the initial 
donor models we use in subsequent calculations. 

We show the distribution in orbital period, P or b, at 
contact (integrated over all M2,i) in Figure [2] with the 
solid line. In this plot, the number, N, in each bin repre- 
sen t the total number of W D channel AM CVn systems in 
the lNelemans et al.l (|200ll ) population synthesis model that 
make contact within each P or b range. For comparison wc 
also show with the dashed line the distribution that results 
from assuming fully degenerate donors. We take P or b at con- 
tact to equal 53.5(7?2,i/O.lP Q ) 3/2 (A'/2,sO.lAf )" 1/2 minutes 
and calculate the fully d egenerate R2,i using a fit to the 
iDelove fc Bildstenl l|2003h zero-temperature He mass-radius 
relation: 



Rl < R2, this tail extends out to 22 minutes, but otherwise 
the contact P or b-distribution is almost unchanged. 

The corresponding ^i-distribution is shown by the his- 
togram in Figure [3] The solid curve in this figure gives the 
cumulative distribution of i/'c.i values. Noteworthy is the fact 
that, unlike the contact P or b-distribution, which is strongly 
biased towards short P or b, the T/>c,;-distribution is approx- 
imately flat between \og(tp Ct i) « 1.4 and 3.6. The differ- 
ence in the two distributions is due to the dependence of 
R2 on tp c ,i'- R2 varies rapidly with ip C: i for ip Ct i ~ 1 — 10, 
but only more slowly for tp c ^ > 100. Thus there is not as 
strong an a priori theoretical preference for large ijj c ^ values 
in the WD channel population as might be expected from 
the contact P or b-distribution. The i/j c ,;-distribution is essen- 
tially unchanged when removing all systems with post-CE 
Rl < R2 from consideration. 



7?2,db = 0.005593 



1.7616 x 10~ 5 
W 2 

+ 0.004643 M 2 



0.008298 In M 2 , (14) 



which is accurate to better than 3% over the range 9 x 
10~ 4 Af Q < M 2 < 0.45M Q . The fully degenerate distribu- 
tion peaks at P or b = 3.0 minutes with a cut-off at 6 min- 
utes. On the other hand, roughly 25% of systems within our 
current calculation make contact at a P or b > 6 minutes and 
5% at P or b > 11 minutes. Additionally, there is a small tail 
of systems making contact out to P or b = 35 minutes (not 
shown in Figure [2j . If we remove all systems with post-CE 



4 THE DONOR'S STRUCTURAL AND 
THERMAL EVOLUTION 

In order to adequately sample the range of initial donor 
properties derived in Sj3]in our subsequent calculations, we 
considered AM CVn evolution starting from 36 different ini- 
tial donor models (solid circles in Fig.[TJ) with M 2 ,; values of 
0.10, 0.15, 0.20, 0.25, 0.30, and 0.325 M Q and log(V> c ,«) val- 
ues of 1.1, 1.5, 2.0, 2.5, 3.0, and 3.5. The coupled binary and 
donor evolution were calculated for two Mi,; for each donor 
model. These Mi,i were chosen so that M to t — M24 + Mi^ 
equalled one of three values — 0.5, 0.825, or 1.325 Mq. 
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Figure 3. The ip c i distribution, integrated over all M 2t i, nor- 
malized to the total number of systems in the population model. 
The solid curve provides the cumulative distribution. Unlike the 
contact P or b distribution in Figure [2] which is strongly peaked 
at shorter periods, the ij) c> i distribution is relatively flat between 
log(i/) c i) « 1.4 and 3.6. 




Log(t/yrs) 

Figure 4. The time evolution of the radius at constant m (black 
lines, the bottom line shows m = O.OOIMq with higher lines 
spaced by A log m = 0.1; the top-most black lines traces -R2), the 
r at which rth = r m (dashed red line, below this line rth > r m ), 
and the lower boundary of the outer convective zone (dashed 
green line). The initial conditions here are M 2 ^i = 0.2 Mq, 
log(V> c ,i) = 2.0, R 2 ,i = 0.0294 R , and M M = 0.3 M . 



To present our results, we first discuss the donor's evo- 
lution in detail for several representative cases. We then ex- 
amine the entire range of our calculations, detailing how and 
in which evolutionary phases different initial conditions af- 
fect the system's parameters. One of our main results is that 
the donor's outer boundary condition is dominated by irra- 
diation from the accretion flow. We also, therefore, present 
additional calculations including the effects of external irra- 
diation and discuss the impact of this additional physics on 
the system's evolution. 



4.1 The AM CVn Donors' Evolutionary Phases 

Figure [3] shows a summary of the donor's evolution in a 
representative calculation. The initial conditions there were 
M 2 ,i = 0.2 M Q , log(Vc,i) = 2.0, R 2 ,i = 0.0294 R e , and 
Mi,« = 0.3 Mq. This model illustrates the range of donors' 
responses to mass loss; for comparison, Figure [S] shows a 
similar summary for a more degenerate, log^c,*) = 3.5, 
donor. In Figure 3] the solid black lines (except for the top- 
most) show the r evolution at fixed m as a functions of time 
(the calculations are started with t = at an a such that 
Rl ~ R2 + 30Hp). The top- most solid line traces out the 
R2 evolution. 

The other lines in this plot indicate various physical 
conditions in the donor. All donor models have, to some 
extent, an outer convective region; the dashed-dotted green 
line indicates the lower radius to which this convective zone 
penetrates. The dashed red line indicates the r at which the 
local thermal time, rth, equals the local mass loss time, r m , 



where these two quantities are defined as 

1* = LT ' ^ , (15) 

r m =£, (16) 

and m! = M—m is the mass exterior to the specified location 
in the donor. Regions below the dashed red line in Figure 
[4] satisfy rth > r m . The evolution in layers where rth S> T m 
is dominated by the nearly adiabatic advection to lower P, 
while layers where rth <C r m are able to adjust their thermal 
structure in response to mass loss almost instantaneously. 

The results shown in Figure 2] are typical of our calcu- 
lations in that the evolution there can be divided into three 
phases. In the first phase, during which M grows from zero 
to its maximum, R 2 decreases towards a minimum value 
in response to mass loss. In Fig. [4] this phase lasts until 
t ~ 2 x 10 5 yrs. The second phase begins as the R 2 evolution 
reverses and the donor begins expanding in response to mass 
loss. This expansion phase corresponds to what is normally 
considered the AM CVn phase of evolution. By this point, 
Tth S» T m throughout the donor so that the donor responds 
adiabatically to the mass loss. Note that the donor enters 
this adiabatic phase well before the R 2 evolution reverses 
(see below) . At some point during the expansion phase, the 
rth = r m line begins moving inward again. This eventually 
leads to the start of the third phase, where the donor is able 
to cool and contract. This is seen at t » 10 9 yrs in Fig. U 
This contraction ends once the donor has shed sufficient en- 
tropy to reach its fully degenerate configuration, even as it 
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Figure 5. Same as Figure [4] but with a more degenerate 
donor; the initial conditions in this system are M 2i i = 0.2 Mq, 
log(V> c ,i) = 3.5, R 2 ,i = 0.0212 R Q , and M lti = 0.3 M© 



continues to cool (evolution beyond t~2 — 3xl0 9 yrs in 
Fig. |U) . 

The first and third phases of the donor's evolution are 
qualitatively new results resulting from our more complete 
treatment of the donor's physics. A more detailed discussion 
of these phases follows. 

4-1.1 The M Turn-on Phase: Evolution to M -Maximum 

Once the WD channel AM CVn donors make contact, our 
calculations show they all begin a phase of radius con- 
traction, a result not anticipated by pri or modelling (e.g., 
iNelemans et alj|200ll ; iDelove et al.1 120051 '). Here, we discuss 
how the donor's structure determines the duration and ex- 
tent of this initial 7? 2 contraction phase. 

The basic picture of the donor's response to mass loss 
can be seen from Figs. [4] and [5] As M 2 decreases, mass el- 
ements move to lower P and p, producing the expansion 
in r at fixed m seen in the evolution of the interior black 
lines. In Appendix [X] we show that this expansion is most 
significant near the surface and that the only contribution 
to R2 contraction comes from the surface term in equation 
(|A2[) . In other words, mass lost from the donor takes with it 
its contribution to R2, tending to produce contraction. The 
overall R2 evolution then depends on whether the under- 
lying layers expand sufficiently to compensate for this lost 
radius contribution. 

Equation (|A5[) shows that this underlying expansion 
depends on the two quantities: \t = (d\nP/d\nT) p and 
V — V, where V' = (d In T/<91n P) m describes the actual 
thermal evolution of a mass element as it is advected to lower 
P. In layers where either of these quantities tends towards 
zero, the p(P) profile remains constant under advection, pro- 
ducing no net change in the relative contribution to R 2 from 
that region. This can occur for strongly degenerate plasmas, 



where \t ~ 0, or when the advected mass element arrives at 
a lower P with the same entropy as the material it replaces. 
Only when the advected mass elements arrive with lower 
entropy (V' > V) is the relative R2 contribution reduced. If 
the latter case dominates in the outer layers, this can lead 
to a net decrease in R 2 . 

Degeneracy effects never dominate in our donors outer 
layers — at least during this contact phase — so the R2 evolu- 
tion depends only on the mode of heat transport in the outer 
layers and on the ordering of rth and r m . When rth <C T m , 
heat transport has sufficient time to redistribute entropy so 
that V' ~ V, and little R2 evolution occurs. Once r m < rth, 
so that V' — > Vad, the R2 evolution begins to depend on 
the background entropy gradient. In convective regions, V 
is essentially equal to V a d, so that adiabatic advection leads 
to minimal local contributions to R2 evolution. In radia- 
tive regions, however, V < V a d; thus only for radiative sur- 
face regions during roughly ad iabatic advectio n will a net 
R2-contraction occur (see, e.g.. lFaulkner|[l976l . for another, 
qualitative, discussion of this). The rate of the R2 contrac- 
tion increases with the entropy of the layer and with the 
steepness of the background entropy profile. 

The structure of our donor's outer layers consists of a 
superficial radiative layer overlying a thin outer-convective 
zone, followed by another radiative region that extends to 
the stellar centre. The thickness of outer-convective region 
increases with donor degeneracy. The entropy profile of the 
inner radiative region is very steep near its outer bound- 
ary; moving inward, this entropy profile tends to flatten out 
rather abruptly. With this general interior structure in mind, 
the i?2 evolution trends seen in Figs. [4] and [5] can be under- 
stood from the above discussion. 

For the log(i/> c ,i) = 2.0 donor (Fig. [3]) , as long as 
Tth -C T m , the R2 evolution is minimal. Once this inequality 
is reversed, the R2 evolution starts to accelerate and a rapid 
contraction ensues due to the steep entropy profile in the 
donor's radiative photosphere. By this time, the underlying 
layers are advected nearly adiabatically, and the entropy gra- 
dient near the surface is continuously decreasing. This leads 
to the decreasing rate of R2 contraction. Finally the entropy 
profile becomes sufficiently shallow that the expansion of the 
underlying layers takes over, and R2 begins increasing. 

The evolution in Fig.[5]is similar, but here the evolution 
is adiabatic (r t h 2> r m ) throughout. The R2 contraction 
rate is slower both because the outer convective zone is 
thicker and the underlying radiative region has both a 
lower entropy and a shallower entropy profile. This also 
results in less overall R2 contraction. Additionally, R2 
decreases significantly only once a sufficient amount of 
mass, 8m, has been lost for the pressure perturbation 
in the outer layers, SP ^ GM2,idm/(iTvR 2 :i ) ~ 9 x 
10 20 dynecm- 2 (5m/M 2ti ) (Af 2 , s /O.lM ) 2 (O.Oli4>/i? 2 ,O 4 , 
to be of order the initial P at the top of the radiative 
region. 

The relative rate of R2 and Rl evolution — characterized 
by £r 2 = (dln7? 2 /dlnM 2 ) and £r l = (d In R L /d In M 2 )— 
determines that of M. In all our evolution models, £r l 3> 
£r 2 initially. However, since 



M 2 



+ q 



GW M2 

£h l rapidly decreases with growing M . Once 



(17) 
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Figure 6. The M(t) evolution for AM CVn systems with M 2l ; = 
0.2 and Mi j = 0.3Mq. Systems with donors having log(i/> Ci i) = 
3.5, 3.0, 2.5, 2.0, 1.5, and 1.1 are shown by the solid, dotted, 
short-dashed, dashed, short-dash dotted, and dash dotted lines, 
respectively. The corresponding R 2 ,i = 0.0212, 0.0225, 0.0248, 
0.0294, 0.0380, 0.0561 Rq, respectively 



these two quantities tend to track each other closely due 
to the sensitive AR dependence of M. Thus, when the £r 2 
evolution is smooth, the M evolution is also. 

This can be seen in Figure [6] which shows the M time 
evolution for a set of models with M 2t i = 0.2, = 0.3Mq 
and differing tp c ,i (indicated by line style with t/> c ,i decreas- 
ing left-to-right). The four lowest ipc^ donors have very thin 
outer convective zones and their corresponding £r 2 evolu- 
tion is smooth. The two highest ip c ,i donors (dotted and 
solid curves) have their R 2 response dominated, at first, by 
their thicker outer convective zones. Once 8m begins prob- 
ing the underlying radiative region, their £r 2 suffer rapid 
changes in slope. This transition from convective to radiative 
dominated R 2 response produces the non-smooth and non- 
monotonic behaviour seen in the dotted and solid curves. 

Figure [6] illustrates other general trends. Initially, Rl > 
R 2 in all cases. While this remains true, the M growth 
is approximately exponential (corresponding to the initial, 
steep increase in M). Lower ip Ci i donors produce slower M 
growth both due to their larger contact a (reducing the 
system's J /J) and larger Hp/R 2 ; the latter tends to be 
the more significant factor. Since the M required to pro- 
duce a > is greater than Mo in all the systems we 
considered, all experience a phase where Rl < R2- Once 
this occurs, the M growth slows as M now depends on 
the donor's non-exponential, sub-photospheric p-profile (see 
equation [13]). This results in the turn-over in M growth seen 



at M 



10" 



-10~ a M© yr" 1 . The final upturn in M be- 



fore each maximum results from a rapid increase in Mo via 
the increasing p of the surface layers. This results from the 
adiabatic advection of lower entropy material to the surface 



and corresponds to the final phase of rapid R 2 contraction 
discussed above. 

The extrema in R 2 , a, and M, which mark the end 
of the turn-on phase, do not occur simultaneously. The 
i?2-minimum happens first and the subsequent expan- 
sion contributes to accelerating the M growth. The a- 
minimum occurs once M reaches a critical value, M cr it = 
— (J/J)gwM 2 /(1 — q) to produce a = 0. The larger J/Jgw 
in systems with more degenerate donors produces higher 
M cr it. Eventually, the a-expansion leads to a decreasing AR, 
and M reaches its maximum. By M-maximum, the donor 
has lost between 5m ~ (0.02 — 0.2)M2,i. Less degenerate 
donors suffer the greater mass loss. Overall the turn-on 
phase lasts between ~ 10 4 and ~ 10 6 yrs. 



4-1.2 The Donor's Cooling Phase 

By the time the system has evolved past M-maximum, the 
donor is evolving adiabatically in response to mass loss. 
Prior work on the evolution of AM CVn systems in this 
phase has neglected the donor's thermal e volution and has 
eithe r relied on predetermined M-R tracks (|Nelemans et al.l 
I2001I ; If armor fc Phinnev|[200l) or has assumed the donors 
conti nue to evolve adiabatically indefinitely (|Delove et al.l 
120051 ). Our current calculations show that at late-times, the 
assumption of adiabatic evolution becomes invalid (as can 
be seen by the evolution of r m = Tth m Fig. 13}. 

Why the donor's adiabatic evolution ends can be under- 
stood by considering how r m and T tri evaluated at m' — M 2 
evolve under mass loss. The quantity 



d In Trr 



din M 2 I m'=M 2 



= 1 



d InM 
d\nM 2 



(18) 



For GW driven, conservative mass transfer, the latter deriva- 
tive can be written as 



d InM 



d In a 



2 - 4 [i R 



(19) 



d\nM 2 d\nM 2 

where the final approximation is good when M 2 <g Mi. In 
the same regime, 

d lnr m I 

=M 2 



dlnM 2 



4 U«2 " o 



(20) 



It is useful to rewrite Tth = J 2 cpTdm' / L as 
c'pT c M 2 /L where c' P = f™ 2 c P Tdm /(T C M 2 ) since both c' P 
and L vary only by a factor of a few during the adiabatic 
phase. In contrast, M 2 and T c vary by roughly one and two 
orders of magnitude, respectively. Holding c' p , L constant, 
d In rth /d In M 2 ~ 1 + d In T c /d In M 2 , and we can use hydro- 
static balance to write 



d lnT c I 
din M 2 I m'=A/ 2 



dlnT c 
d\nP c 



d lnP c 
dlnMo 



= V ad , c (2-4£fl 2 ), 
(21) 

where V a d, c characterizes the T c evolution during the adia- 
batic phase. 

Thus once M 2 <C Mi in the adiabatic phase, r t h/T m oc 
M 2 where 

10 

3 

Typically, £ R2 w -0.3—0.2 and V ad , c ~ 0.36-0.39, so that 
a ~ —5.7-— 5.2 and leading a rapid evolution of r t h/T m as 



4£a 2 



+ V ad , c (4£fl 2 - 2) 



(22) 
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Figure 7. The evolution of T t h and r m (at m' = M 2 ) as a function 
of M 2 for the set of AM CVn models with M 2 ,j = 0.2M Q and 
Mi j = O.3M0 shown in Fig[6] The various lines styles have the 
same meaning as in that figure. 

M 2 is reduced. In the early portion of the adiabatic phase, 
where M 2 ~ 0.1M Q , r th /r m ~ 10 3 - 10 5 and this scaling 
then gives T t h/V m ~ 1 when M 2 ~ 0.01 — O.O3M0. 

We compare this prediction to our numerical calcu- 
lations in Fig. [7] where we show the computed evolution 
of both time-scales for our models with M 2; i = Q.2Mq, 
M hi = O.3M . The ratio r th /r m w 1 at M 2 = 0.01 - 0.02 
in all cases, in excellent agreement with the scaling rela- 
tion. Beyond this point, the donor is able to shed entropy, 
allowing less degenerate donors to contract towards a fully 
degenerate configuration. 

4.2 Impact of Initial Conditions on Donor 
Evolution 

We now explore how the donor and binary evolution varies 
with initial conditions. In Fig. [5] we show the i? 2 (M 2 ) evo- 
lution for several sets of calculations. In this figure, differ- 
ent colours correspond to different initial donor degeneracy, 
while different line styles indicate different Mtot • The evolu- 
tion along each track is from right to left, with the donors 
first evolving steeply downward in i? 2 . The amount R2 de- 
creases prior to the i? 2 -minimum depends on ip c ,i- Less de- 
generate donors contract to a greater extent during the turn- 
on phase as discussed in £14.1.11 

In panel (a) of Fig. [5] the R2 evolution for a set of sys- 
tems with Mtot = 0.825Afo are displayed. During the turn- 
on phase and early expansion phase, i? 2 (M 2 ) depends on 
both ip Ct i and M 2i i. At fixed i\i Ct i different M 2 ,i have differ- 
ing s(m) profiles, and this is reflected in the R2 evolution. 
The entropy differences are more significant at larger m, pro- 
ducing the tendency towards convergence as M 2 is reduced. 
Near the local _R 2 -maximum (at M 2 « 0.01 — 0.03), track 
convergence is furthered by donor cooling, which erases any 




log(M 2 /M ) 

Figure 8. The dependence of the i? 2 -evolution on initial binary 
parameters. In both panels, colours indicate the donor's initial 
log(-0 Ci i): 1.1 (yellow), 1.5 (red), 2.0 (green), and 3.0 (blue). Line 
style indicates the binary's Mtot: 0.500 (dashed), 0.825 (solid), 
and 1.325 Mq (dot short-dashed). Panel (a) compares the _R 2 - 
M 2 evolution for systems with Mtot = O.825M0 to examine the 
impact of initial degeneracy and donor mass. We show evolution 
for systems with M 2j i = 0.15, and 0.3 Mq and the initial de- 
generacy listed above. Panel (b) examines the impact of varying 
Mtot on the R2 evolution. We show comparisons at two i\i c ,i to 
illustrate how the magnitude of Mtot effects are reduced at higher 
donor degeneracy. 

remaining initial s-profile information. By this point, the 
donors cool and contract along tracks parametrized by ip c ,i- 
Cooling continues and R2 contraction slows as the donors 
become increasingly degenerate and finally reach the fully 
degenerate M-R relation. 

In panel (b) of Fig. [8] we consider the dependence of 
the R2 evolution on Mtot- During the turn-on phase, the 
donor evolution is essentially independent of M to t ■ Only for 
the very lowest degeneracy donors is this not the case. These 
donors are able to cool somewhat during the turn-on phase 
and systems with lower M to t produce slower M growth, al- 
lowing greater cooling before mass loss becomes adiabatic. 
The larger M in higher M to t systems also extends the du- 
ration of the adiabatic expansion phase to lower M 2 and 
larger _R 2 . The R2 path a system follows during its later 
contraction phase is thus parametrized by i[> c ,i and M to t- 

In Fig. [9] the evolution of the donor's L is shown as a 
function of P or b , with evolution proceeding downward along 
each track. The L-evolution reflects the donor's evolution- 
ary stages. The significant decrease during the turn-on phase 
(while Porb is decreasing) is due primarily to the removal of 
mass from the steep s-profile region (see Appendix|B]). Once 
the steep s-gradient material has been removed, the L evo- 
lution slows along with the R2 contraction. In the adiabatic 
expansion phase, L roughly plateaus at a level set by the core 
entropy profile. That is, L during this phase is parametrized 
by ip c ,i and Porb- Variations in the s profiles due to M 2j i 
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Figure 9. The dependence of the donor's L-evolution on initial 
binary parameters. The line styles and colours have the same 
meanings as in Fig. [8] As in that figure, panel (a) shows the 
effects of varying M^^ and tp C: i at fixed Mtot, while panel (b) 
examines the effects of varying A'Lot, M2 : i at fixed if) c ^i. Here M2,i 
has values of 0.f5, and 0.3 Mq. Evolution proceeds downward 
along each track and M^^i can be inferred from the initial P or b 
within each set of constant tp c ,i tracks, with Mi,i decreasing with 
increasing P or b i- 
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Figure 10. The dependence of the M-P or b relation on initial 
binary parameters. The line styles and colours have the same 
meanings as in Fig. [8] For clarity, we display the evolution of 
two (Mi,*, M 2 ,i) pairs: (0.350, 0.150), and (1.025, 0.300), where 
values are in Mq and values of log(i/> Cj i) of 3.0 and 1.1. The figure 
focuses on the M evolution post M-maximum. Evolution is first 
upward along the vertical sections of each track, then outward in 
P or b as M decreases. The range of M at fixed P or b displayed here 
represents the range available to WD channel AM CVn binaries 
given our assumptions in determining the donor's initial entropy. 



differences are reflected in slight differences between tracks 
sharing ip ,i- During this phase, L ~10 _6 -10~ 4 Lq, while 
T e g also plateaus at values «1000-1800 K. The slow evo- 
lution in L ends once donor cooling sets in, producing the 
rapid decline starting between P or b ~40-50 minutes. Panel 
(b) shows that a larger M to t extends the P or b at which donor 
cooling start (via an increased M); M to t thus parametrizes 
the L evolution during the cooling phase and beyond, with 
higher Mtot producing higher L at fixed P or b- 

The differences in ^2(^2) evolution with M to t and ip Cl i 
are reflected in the binary's M(P or b) evolution as shown in 
Fig. 1101 For clarity, the evolution is shown for only four sys- 
tems; these systems do, however, show the range of phase- 
space covered by the overall population in this study. We 
show models with \og(ip Ct i) —3.0 and 1.1, the latter provid- 
ing a reasonable lower limit to the degeneracy expected in 
this population. 

Although not shown by this plot, the M-P or b evolution 
during turn-on is determined primarily by Mi,t and ip c ,i\ 
Mtot influences only the M maximum at large ip c ,i- Once 
the systems have evolved into the expansion phase, the P or b- 
M evolution depends only on M to t and ipc,i, since donors 
with the same %j} c ^ eventually follow the same R^i^h) evo- 
lution. For GW-driven j, at fixed P orb , M oc Mf^M? (for 
M2 <C M2) so that the variations in the donors' R2 evo- 
lution produce a larger change in M than do M to t vari- 
ations. This trend is shown in Fig. 1101 where the family 
of curves with \og(tp Ct i) — 1.1 lie significantly above those 
with \og(ip Ct i) = 3.0. Within each set of curves at fixed ip c ,i, 



the spread in M is caused by M to t differences. Once donor 
cooling sets in, all tracks collapse toward the fully degener- 
ate ones. Systems with lower ip Ct i and larger M to t reach the 
fully degenerate tracks at longer P or b- After this point, the 
Af-Porb relation is parametrized exclusively by M to t- 

4.3 Evolution with Irradiative Feedback 

4-3.1 Comparison of the Irradiative and Donor's Intrinsic 
Fluxes 

The donor's thermal evolution does not occur in isolation, 
but in the radiation bath provided by the flux from the 
accretor and disk. As the compressional heating luminosity 
generated in the accretor is always m uch less than the accre- 
tion luminosity (jBildsten et al.ll2006l ), the external flux seen 
by the donor is dominated by the accretion light. To deter- 
mine how important this external flux may be to the donor's 
evolution, we calculate the accretion luminosity, Lace, using 

Lace = M((j> L1 -(pRl), (23) 

where </>li, 4>ri are the gravitational potential at 
the in ner Lagrange point and the accretor, respec- 
tively (|Han fc Webbinkl 1 19991 ), and then compare Ti rr — 
(Lacc/ 47rq-a 2 ) 1 ' 4 to the donor 's T e s. We use equations (14)- 
(16) of lHan fc Webbinkl (|l999T ) to calculate 4>l\, <f>m (noting 
these authors interchange our definitions of Mi and M2). 

The evolution of T e g/Ti Tr versus T e g for the evolution 
tracks from Figure [10] is shown in Figure 1111 During the 



© 0000 RAS, MNRAS 000, 000-000 



The Thermal Evolution of the Donors in AM CVn Binaries 11 



0.5 



-0.5 



M 
O 



1.5 



' \ 



I 1 I 

/ 

I 

J / 

/ 1 



\ / 



_^ I I 

3.5 

log(T eff /K) 



Figure 11. A comparison between the T e ff of non-irradiated 
donors and the Ti rr produced by the system's accretion flow. We 
show the same set of models as in Figure 1101 Line styles and 
colours have the same meaning here as in that Figure. For the 
entire range of donors, irradiation dominates the donor's outer 
boundary condition by the time M has grown to its maximum 
value. 



turn-on phase, Tj rr increases rapidly, and by M-maximum, 
Ti rr has grown larger than T e ff by a up to a factor of 50. In all 
but the lowest Mtot and ip Ct i cases, Ti rr remains greater than 
T e ff after this point. In most cases, T[ TI strongly dominates 
T e ff. We consider now how this fact impacts the donor's and 
binary's evolution. 



4-3.2 Our Irradiation Modelling 

The general effect of external irradiation i s to increas e 
the temperature of the donor's atmosphere (|Milnd fl92rj ). 
whic h will tends to low e r L and slow its cooling (see, 
e.g., iBurrows et~ai1 120031 ; iBaraffe et ail 120031 , in the con- 
text of irradiated planets). The amount of heating de- 
pends on many factors: the irradiating flux's intensity, 
spectrum, anisotropy, and the opacity sources in the 
donor's atmosphere that determine its a l bedo (see, e.g., 
Vaz fc Nordlundlll985l;lBarman et all 1200 ll; ISudarskv et all 
20031 ; iBarman et al.ll2004l ; lBurkert et alj|2005h . The detailed 
modelling of all these effects involves a multidimensional, 
non-grey, radiation-hydrodynamics problem and is beyond 
the scope of this paper. Thus, to model irradiation's impact, 
we continue to assume grey opacities and alter our temper- 
ature outer boundary matching condition to T e = T p hot, 
where T p hot is defined by 



rpi rp4 . Thr 1 ( L . Lace 



(24) 



where we assume a point-source geometry for L acc and com- 
plete red istribution of the ir radiating flux around the donor's 
surface (jRitter et al.l l2000h . We also define the quantity 



L BU r! = 47T(ji?2Tp hot , which gives the total (intrinsic plus 
thermalized irradiative) luminosity off the donor's surface. 
The factor 77 is a dimensionless efficiency parameter giving 
the fraction of L acc that is thermalized in the donor's pho- 
tosphere. Equation (|24|l makes it clear that 77 parametrizes 
(within a grey, ID model) our entire ignorance associated 
with the above uncertainties (including uncertainties in the 
assumed geometry). 

For our calculations, we choose a fixed 77 < 1.0 and 
calculate L acc from the system's secular M. However, AM 
CVn binaries experience a phase in which instabilities in 
the He accretion disk p roduce cyclical variations in Mi 
|Tsugawa fc Osakl ll997T l. While these outbursts produce 
brightness v ariations < 4 magnitudes, the outburst perio d 
is ~ 5 days (|Wood et alj|l987t IPatterson et al.ll 19971 . l200Ch . 
This corresponds roughly to a r t h at ml < 10 _10 M2, so al- 
most the entire donor is only aware of the time-averaged 
iacc, making our use of the secular M reasonable. 



4-3.3 Impact of Irradiation on Donor & Binary Evolution 

A star's outer boundary condition only has a significant im- 
pact on it s structure when the star is ne arly fully convective 
(see, e.g.. lKippenhahn fc Weigertj|l99d Sec. 10.3). Further, 
during the adiabatic phase of AM CVn evolution, mass loss, 
not thermal processes, dominate the donor's evolution. Thus 
we can expect that irradiation will affect the donor's evo- 
lution most during the cooling phase, by which time the 
donor's are fully convective and thermal processes dominate 
their evolution. This expectation is borne out by our numer- 
ical calculations. 

To illustrate this, we show in Figure [12] a comparison 
of T(P) profiles between an irradiated (77 = 0.5, red lines) 
and non-irradiated (77 — 0.0, blue lines ) donor at several 
specified values of M2 along respective binary evolution cal- 
culations. Apart from the 77 differences, the initial condi- 
tions for the two calculations are the same. The profiles 
for M 2 = 0.1999 (solid lines), 0.1998 (short-dashed) and 
0.1904Mq (dashed) all occur during the turn-on and adi- 
abatic phase. By M2 = 0.015Mq (short-dash dotted line), 
the non-irradiated donor is in its cooling phase. 

Throughout, Tin- dominates T p hot, so that the 77 = 0.5 
donor's T p hot-evolution tracks that of M. By increasing 
T P hot , a constant Tj T I will t end to reduce L iBaraffe et al.l 
120031 ; IBurrows et"al] 120031 ; I Arras fc Bildstenl I2006T ). Here, 
there is at least one other effect important to the L evo- 
lution. Irradiated envelopes tend to have more extensive ra- 
diative regions and steeper entropy profiles. During rapid 
mass loss, this increases the L decrement a mass element 
experiences as it is advected outwards. 

Both effects contribute to lowering the donor's L in the 
irradiated model. During phases of rapid mass loss, the lat- 
ter effect even leads to a net L < (most obviously seen 
in the inverted T(P) profile of the M 2 = O.19O4M case, 
but also present at M 2 = O.1998Af ). A net L < is pro- 
duced when the flux cost required to advect material up the 
steep entropy gradient cannot be provided by donor's intrin- 
sic flux; the deficit is made up by absorption of irradiating 
flux and an inverted T(P) profile results. Although some of 
the irradiating energy is absorbed below the photosphere, 
never more than 10 _6 M2 of the donor is involved and the 
inward directed flux is always < O.OI77 o T-^. The persistence 
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Figure 12. The T(P) profiles for anij = 0.5 irradiated (red lines) 
and non-irradiated (blue lines) donor at points in their evolution 
when M 2 = 0.1999 (solid lines), 0.1998 (short-dashed), 0.1904 
(dashed), and O.O15M0 (short-dash dotted); dotted portions of 
the lines indicate convective regions. The initial conditions for 
both sets of calculations were M 2 ^ = 0.2, Mi i = O.625M0, and 
R 2i i = 0.0380-Rq. Only once the non-irradiated donor is mainly 
convective and has begun its cooling phase (M 2 = 0.015Mq) 
do significant interior structural differences with the irradiated 
model appear. 



of the L < condition depends on M and rj, with L recover- 
ing to positive values once M decreases sufficiently. For the 
case shown in Fig [T2] L < until M 2 ~ 0.02M Q , beyond 
the point the non-irradiated donor started its cooling phase. 

Figure [12] further illustrates our general results that dif- 
ferences between the deep interior structures of irradiated 
and non-irradiated donors only extend as far as the base 
of the convective layer, roughly speaking. Thus it is only 
once the donors become nearly fully convective do signif- 
icant structural differences occur; such differences become 
even more apparent once the non-irradiated donor begins 
its cooling (e.g., the M 2 = 0.015M© profiles in Fig. [HJ. 

By decreasing L, the net effect of irradiation is to ex- 
tend the adiabatic phase of evolution to lower M 2 and to 
slow the donor's cooling afterwards. Irradiated donors con- 
tinue expanding to longer P or b and contract more slowly 
once cooling begins. This alters the M evolution, as illus- 
trated in Figure [T3] The solid lines show systems' evolu- 
tion with non-irradiated donors, the short-dash dotted and 
dashed lines show rj = 0.1 and rj = 0.5 tracks, respectively. 
It can be seen that differences in M evolution between rj val- 
ues appear even before the cooling phase for non-irradiated 
case. This reflects the outer boundary condition's growing 
importance in r\ = 0.0 donors, which have much deeper con- 
vective regions by this point in the evolution. As rj is in- 
creased, the P or b where the M-decline occurs increases; this 
effect is more pronounced at lower ip c ,i- Before contracting 
to the fully degenerate track, irradiated donors converge to 
an 77-dependent, intermediate cooling track. In Fig. 1131 this 
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Figure 13. The impact of irradiation on the binary's M evolution 
for the initial conditions M2 j = 0.2 and = O.625M0. Line 

colour indicates log(i/> Ci i): 1.1 (red), 2.0 (green), and 3.0 (blue); 
line style indicates rj: 0.0 (solid), 0.1 (short-dash dotted), and 0.5 
(dashed). Only the non-irradiated log(f/) c i) = 3.0 case is shown 
due to numerical difficulties in converging the irradiated models 
during the turn-on phase. However, from other irradiated models 
with log(i/! Ci i) = 3.0 the irradiated log(i/> Cj i) = 3.0 tracks will 
not differ substantially from the r\ = 0.0 track shown. For lower 
ip c i, irradiation extends the adiabatic phase and slows the donor's 
cooling, elevating M during these phases. 
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Figure 14. Irradiation's impact on the donor's L SUI f in the limit 
of zero Bond albedo for systems sharing the initial conditions of 
M 2 ,i = 0.2, Mi,; = O.3M , and log(»/>c,i) = 2.0. Line style has 
the same meaning as in Figure [13] In the r\ 7^ cases, irradiation 
substantially elevates L sur f. 
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occurs at P or b ~ 49(54) minutes for the 77 = 0.1(0.5) tracks. 
By this point, the irradiated donors are fully convective, so 
that R2 depends on s c , M2 , and the strength of irradia- 
tion (| Arras fc Bildsten|[200r3) . As irradiation is a function of 
M-i{R2) through M, the evolution of fully-convective donors 
within our model are parametrized exclusively by M to t and 
77, as evidenced by this intermediate convergence. 

The evolution of L sur f for a set of donors with differing 77 
at fixed Mi,i, Ma.i, and ip c< i are shown in Figure [141 Within 
our grey-atmosphere modelling, Fig. Q3] can be interpreted 
as providing the donor's surface luminosity in the limit of 
zero Bond Albedo. Since Tj rr dominates T p hot over most of 
the evolution, the L sur f evolution mirrors that of M. Dur- 
ing the turn-on phase, L SUI f decreases initially in all donors. 
This decrease is reversed for 77 7^ by the M growth. By M 
maximum, L suri ~ 10 3 -10 5 higher in irradiated models com- 
pared to non-irradiated donors. During the adiabatic phase, 
Z/ SU rf decreases with M , while the non-irradiated donor's L 
plateaus. In the cooling phase L BUT { converges towards tracks 
parametrized by 77, but along much shallower slopes than 
non-irradiated donors. Even so, by P or b ~ 60 minutes, even 
77 = 0.5 donors are still rather dim, with L sur f < 10 _6 L©. 



5 DISCUSSION AND APPLICATIONS 

5.1 Comparison to Prior AM CVn Donor Models 

We start our discussion by comparing the M(P or b) evolu- 
tion produced in our current donor models to that pro- 
duced by several prior donor models in Figure 1151 The 
lower panel focuses on the c omparison between our c urrent 



models (solid lines) and the iDeloye fc Bildsten" d2003l) isen - 

](|20 

tion of these mode ls in both IDeloye fc Bildsten] (|2003l ) and 



iye 

tropic models (dashed lines). The lDelove fc Bildsten! (|2003l ) 
models assume the donors ar e fully convective and app lica- 



iDelove et al] (|2005l ) also assume adiabatic donor evolution. 
Our more complete, current modelling shows that these He 
donors are not fully convective over much of the AM CVn 
evolution phase and that adiabatic evolution only occurs out 
to P or b ~ 40 — 55 minutes. 

The impact of these differences is apparent between the 
M evolution shown by the dashed and solid lines in Fig. 
1151 The initial conditions (i.e. Mi, M2, R2) for each dashed 
line evolution equals the set of these values along the corre- 
sponding solid line at the point of intersection. From there, 
the isentropic donors produce M evolution that increasingly 
diverges upwards from the solid-line tracks, even during the 
adiabatic evolution phase. This results from the realistic 
models being radiative throughout much of their core. Thus, 
compared to isentropic donors, the realistic models have a 
greater £r 2 < 0, even during adiabatic evolution ( i]4.1.1l and 
AppendixfX]| . The differences between these evolution tracks 
become even more dramatic once the realistic donors be- 
gin cooling. Thus, modelling ultracompact binary evolution 
using isentropic donors will overestimate the R2 expansion 
rate, the M-P or b relations, and P or b(i)- 

This has implications for how M measu rements con- 
strain an AM CVn system's formation channel. IDeloye et al.l 
(|2005l ) showed that assuming isentropic donors leads to a 
significant overlap in the M-P or b plane between hot WD 
channel systems and He-star channel systems. The overlap 
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Figure 15. Comparison of AM CVn system's outward Af-P G rb 
evolution produced by different donor models. The lower panel 
compares evolution with our c urren t models (solid lines) to that 
with the IDeloye fe Bildsten! J2003T) isentropic models (dashed 
lines). The two solid lines had M%,i = 0.2Mq and log(i/> c = 1.1 
(upper) and 3.5 (lower); in all four systems Mtot = 0.5Mq. The 
initial conditions [M\, M2, -R2) for each dashed line evolution 
equalled the set of these values on the corresponding solid line 
at the point of intersection. The upper panel compares evolution 
with our current models (solid lines, lower same as in lower panel, 
upper with Aftot = 1.325, M 2 ,j = 0.3M o , log(i/> c ,i) = 1-1) with 
the ev olution produced by th e semi-degenerate donor M-R rela- 
tion of lNelemans et all 1120011 ) with Mi,, = 0.4 and M 2 ,i = O.3M 
(red line). 



between these channels is much reduced by with our cur- 
rent donor models, as illustrated in the upper panel of Fig. 
1151 There the black lines show the maximum range in the 
Af-P or b relation for WD channel AM CVn systems produced 
with our current unirra diated models. The red line shows the 
iNelemans et al.l (|200ll ) fit to a semi-degenerate donor M{R) 
evolution with Mtot = 0.7Mq, and provides a rough lower 
limit to He-star channel systems' M(P orD ) evolution. Only 
the most massive WD channel systems have any overlap 
with He-star channel systems. Thus, given this current the- 
ory, a determination of a secular M significantly above the 
upper black Fig. [l5]would indicate a system formed through 
the He-star channel. This statement has one caveat: if WD 
channel systems have a maximum entropy greater than indi- 
cated by our determinations in ij3]or if the donors are heated 
significantly (e.g., by tidal mechanisms) earlier in the AM 
CVn phase, then the overlap with the He-star channel sys- 
tems could be increased. 

We point out that mo delling He-star channel systems 
with INelemans et al.l (|200lf ) semi-degenerate M(R) relation 
beyond P or b > 45 is problematic. This is because these 
donors will also begin cooling and contracting by these P or b, 
similar to the WD channel donors. Thus, in reality, the red 
line in Fig [15] should begin a down turn somewhere in the 
vicinity of the Mtot = 1.325M© track. In fact, the original 
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calculation used by_ Nelemans et al.l (12001 ) to determine this 
fit (model 1.1 of lTutukov fc Fedorovall 19891 . see their Fig. 2) 
shows the start of this down turn at P or b ~ 40 minutes. 
Thus, beyond the P or b at which donors begin their cooling, 
M measurements alone will not distinguish between forma- 
tion channels. 



5.2 The Orbital Period Distribution of 
WD-Channel AM CVn Binaries 

While the changes to the adiabatic phase M 2 (R 2 ) evolu- 
tion will quantitatively alter our expectations for the AM 
CVn population's P or b-distribution, the detailed evolution 
to Porb-minimum and the occurrence of donor cooling will 
produce qualitatively new features in the population's Porb- 
distribution. We discuss these features here. 

A key factor in how an individual system contributes 
to the P or b-distribution is the time-derivative, P or b, of its 
Porb-evolution. In a steady-state, continuity requires that 
the number density of syste ms at some P or b, np nrU , s cale as 
n P„ b K IV-Porb| (see, e.g.. iDelove fc Bildstenl 120031 ) . Here 
we will use this scaling relation to display steady-state dis- 
tributions of systems sharing initial data. This provides a 
straightforward means for displaying how differing initial 
conditions influence the relative contribution systems make 
to the overall P or b-distribution. We calculate P or b along an 
evolutionary track via 



Porb 
Porb 



GW 



(25) 



We show the evolution of P or b about the P rb-minimum 
for a representative set of systems in Figure [16] These sys- 
tems have Mi,i = 0.575, M 2t i = 0.25Mq, with colours indi- 
cating different tp c ^. The short-dash dotted segments indi- 
cate Porb < 0, solid segments P or b > 0. The black line shows 
the inward evolution due only to the ( J / J)gw term in equa- 
tion (|25p for comparison. As discussed in ^4. 1.11 the evolu- 
tion here is most sensitive to M 2t i and tp c ,i - The P or b evolu- 
tion reflects this, with ip Ct i (as illustrated in this figure) and 
M 2 ,i both affecting the the value of the P or b minimum and 
how Porb diverges from the GW-only evolution once contact 
occurs. The generic features of our this P or b evolution are 
the strong spikes at P or b-niinimum as P or b evolves through 
zero and the existence of P or b minimum and maximum that 
occur before and after P or b-minimum, respectively. 

The relative magnitude of P or b before and after the 
Porb -minimum — which determines the relative number of 
systems evolving inward vs. outward in a steady-state — 
depends on M 2ti and ip Cii . The P or b < just before Porb- 
minimum is approximated by the GW-only rate, which 
scales as (J/J)gw- By P or b-maximum, M has achieved its 
secular rate: 



Ma 
M 2 



{J / J)gw 
5/6 + W2-q' 



(26) 



producing a P or b during expansion: 



Porh 



GW 



~ 1/3 



Ch 2 + 5/3 - 2q 



(27) 

If M 2 does not change appreciably from P rb-minimum 
to Porb-maximum, the term in square brackets, denoted here 
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Figure 16. The evolution of P Q rb about the P rb- m i n i mum f° r 
representative systems with M 2 i = 0.25 and Mi 4 = 0.575Mq. 
Line colour indicates log(^) c j): 1.1 (yellow), 1.5 (red), and 3.0 
(blue). Short-dash dotted line segments indicate P < 0, solid lines 
indicate P > 0. The black line shows the corresponding inward 
Porb evolution produced by GW emission alone. 



by Pp, estimates the relative magnitudes of P or b before and 
after P or b-minimum. Note that pp depends on ip Ci i through 
£r 2 . During the early expansion phase, £u 2 evolves from to 
~ —0.3 to —0.05. With < q < 2/3 (the q range considered 
here), Pp can vary considerably. Typically, — 1 < /3 p < 0, as 
seen in Fig. \W\ However, for larger q and tp c ,i, Pp can be 
much less than -1; we see this behaviour in our calculations 
with Mi,, = 0.3, M 2 ,i = 0.2M Q for log(V> c ,i) > 3.0. Lower 
ip Ci i donors have smaller relative outward-to-inward |P or b|, 
for two reasons. One, lower ip Ct i leads to larger minimum 
£h 2 , producing slower expansion and Pp values closer to 
zero. Two, M 2 during the P or b-transition is not fixed, with 
M 2 losses greater for lower ip c ,i- This results in a lower J 
loss rate post P rb-maximum, further contributing to smaller 
Porb > 0. 

How the evolution near P or b-minimum translates into 
the steady-state P or b-distributions is shown in Figure 1171 
There we display histograms of relative system numbers 
along evolutionary tracks with Mi,j = 0.575, M 2l i = 0.25. 
The number in each bin, N, is calculated by integrating 
np oTb dP OI h over each bin. The overall normalization is ar- 
bitrary but fixed across tracks, so that the histograms ac- 
curately reflect relative P or b rates. The lower panel displays 
separately the contributions from the P or b < (short-dash 
dotted lines) and P or b > (solid lines) segments of the evo- 
lution. The black line in the lower-panel shows the corre- 
sponding result for a GW-only driven in-spiral. The upper 
panel shows the sum of contributions from both inward and 
outward evolving systems. 

The general trends in the lower panel are the slight de- 
viation from GW-only evolution from contact inward, fol- 
lowed by a peak as P or b — > at P or b-minimum. The out- 
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Figure 17. Histograms of relative numbers of AM CVn sys- 
tems in a steady-state along single evolution tracks. All tracks 
have Mi^i = 0.575 and M2 i = 0.25Mq with colour indicating 
log(Vc,i): 1.1 (yellow), 1.5 (red), 2.0 (green), 3.0 (blue), and 3.5 
(cyan). The lower panel separates systems by P or b sign: the short- 
dash dotted histograms counts inward moving systems, the solid 
histograms outward moving systems. The upper panel displays 
the sum of both inward and outward moving systems. The black 
short-dash dotted line shows the distribution for the evolution 
produced by GW emission alone (i.e. M = 0). The bins have 
width of AP or b = 6 s. The overall normalization is arbitrary, 
but the relative normalization accurately reflects the relative P or b 
rates. 



wardly evolving systems likewise show peaks just post P or b- 
minimum. For larger P or b, lower-i/) Cj i tracks lead to fewer 
systems at a given P or b sin ce hotter donors ha ve larger M2, 
producing larger ( J/J)gw (|DeIove et alj|2005l . see also Fig. 
I16[) . The sharp steps in the upper-panel histograms result 
from starting our calculations at the point of contact (i.e. 
pre-contact evolution is not included), so the location and 
size of each step is somewhat artificial. However, since they 
result from having a definite starting point for the GW- 
driven in-spiral, as is provided by the CE-event forming 
these systems, there is a physical basis for expecting their 
existence. 

The most striking feature of these histograms are the 
strong peaks at each system's P or b-minimum. The impact of 
these features on the integrated AM CVn P or b-distribution 
below P OI .b ~ 15 minutes will depend on the distribution 
of in itial conditions, th e survival of systems at contact (see, 
e.g., iMarsh et alj|2003 ). and how He-star channel systems 
contribute in this P or b range. Determining how the inte- 
grated Porb-distribution depends on initial conditions and 
the physics determining the outcomes at contact is the sub- 
ject of current work. This will be most relevant to future 
space-based GW- interferometers, such as LISA, which will 
provide a rather complete census of the galactic AM CVn 
population at these P or b (|Nelemans et al.l |200ll ) and offer 
direct observational tests of these predictions. 



In iDelove fc Taaml (|2006h . we considered the relevance 
of the Porb < phase to the short period X-ray variables, 
RX J0806+1527 (32 1 s; iBeuermann et all Il999| ) and RX 
J1914+2456 (569 s; lHaberl fc Motchl 1 19951 ). There is still 
some questio n as to the nature o f these two sources (see 
discussion in lDelove fc Ta am 2006], and references therein). 
If bo th system's periods ar e orbital and their measured 
Porh JStrohmaverll200 J . 120051 ) secular, then lDelove fc Taaml 
(|200€ ) showed that both are consistent with being members 
of the AM CVn population. If this is the case, there is the 
question of whether observing either system in its present 
state is an extremely unlikely event. 

A detailed answer to this question is beyond the scope 
of this paper. However, a related question can be posed: 
what is the relative likelihood of detecting each system in its 
Porb < versus Porb > phase? Figure [l8]provides a partial 
answer. The lower panel shows the P or b evolution for two 
sets of AM CVn systems, each of which bracket the evolution 
of systems consis tent with either RX J0806+1527 or RX 
J1914+2456 (see IDelove fc Taaml liblirj ) . The upper panel 
shows, for these evolution tracks, the ratio of the number of 
systems with P or b > to those with P or b < in steady-state. 
This ratio is only ~ 1.3 — 2.3 for periods below each system's 
measure P or b- Thus, there is not a strong a priori bias to 
detecting such systems in their P or b > phase, although the 
actual relative detection likelihood between the two phases 
is likely influenced strongly by selection effects. 

We now turn to how donor cooling influences the Porb- 
distribution. By the time P or b has passed its maximum, it 
is evolving at a rate given by equation (|27|l . and is essen- 
tially determined by M (which is a function of Mtot and 
ip c ,i)- This is seen in the lower panel of Figure 1191 where 
we display P or b versus P rb for two different M to t at two 
different ip c ,i- Before systems begin their cooling phase, the 
ordering of tracks is for lower tp Ct i donors to produce higher 
Porb- After the cooling phase ends, tracks are distinguished 
only by Mtot- During cooling, the donor's contraction stalls 
the Porb evolution, producing a reversal in the ordering of 
Porb with tp Ct i and a distinctive peak in the steady-state P or b 
distribution along each track (upper panel of Fig. I19|l . 

The location and magnitude of these peaks are deter- 
mined by Mtot and ?/> c ,;. The system's M to t determines the 
Porb at which cooling becomes important, with larger Mtot 
moving the peak's centre to longer P or b- The donor's il> c ^ de- 
termines the degree of donor contraction; hotter donors lead 
to slower P or b during the cooling phase and larger peaks. In 
slowing donor cooling, irradiation also acts to shift the peak 
centres to longer P or b- The distribution of AM CVns above 
Porb ~ 40 minutes will provide an integrated diagnostic of 
the distribution of 4>a, Mtot, and rj In this P or b range, AM 
CVn binaries will not be individually resolvable GW-sources 
due to the galactic fore ground of detached WD- WD binaries 
l|Nelemans et al.ll200ll ). Thus, whether the system distribu- 
tion in this Porb range can be a practical diagnostic tool must 
await observations progress in the optical/IR wave bands. 



5.3 The Eclipsing AM CVn SDSS J0926+3624: 
Evidence for Non-Zero Entropy Donors from 
the WD Channel 

Although we reserve detailed applications of the models pre- 
sented here to specific AM CVn systems for a later compan- 
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Figure 18. The evolution rates for AM CVn system models 
consistent with measured properties of RX J0806+1527 (lines at 
P or b ~ 5 minutes) and RX J1914+2456 (lines at P or b ~ 9 min- 
utes). The set of two lines shown for each system indicates the 
approximate range of variation allowed by observations and our 
modelling. Different colours are meant only to guide the eye in 
matching tracks within and between the two panels. The lower 
panel shows the P or b evolution for the selected tracks, with solid 
lines indicating P orD > and short-dash dotted lines P or b < 0. 
The upper panel shows the ratio of of number of systems with 
Porb > to those with P or b < assuming steady-state. The dot- 
ted lines indicates the measured P or b of each system, while the di- 
amonds show their measured P or b. In making these comparisons, 
we explicitly assume the measured P or b reflects the system's sec- 
ular P or b. 



ion paper, we shall discuss the recently di scovered eclips- 
ing A M CVn binary, SDSS J0926+3624 (|Anderson et al.1 
120051 ) as a system for which the entropy of the donor can be 
probed. This is especially important as the distribution of 
donor entropy provides (i) a potential diagnostic of AM CVn 
formation channels and (ii) constrains the stell ar binary evo- 
lution determining AM CVn initial conditions (|Delove et al.l 
l2005l ). lBildsten" et al.l (|2006T ) have presented evidence based 
on accretor properties that the two known long-period AM 
CVn binaries, CP Com (P rb = 46.6 minutes) and CE-315 
(Porb = 65.1 minutes) both harbour a relatively hot donor. 
Additionally, if early-contact AM CVn binaries are indeed 
the correct model for the systems RX J0806+1527 and RX 
J1914+2456, then this would provide further evidence for 
hot donors: the bracketing models shown in Fig. 1181 have 
a range of log(V> c ,i) ~ 1.3 - 1.8 for RX J0806+1527 and 
« 1.2 - 1.5 for RX J1914+2456. SDSS J0926+3624, how- 
ever, has finally provided direct evidence for hot donors in 
AM CVn binaries. 

This evidence comes in the form of Mi and M 2 de- 
terminations mad e via modelling of eclipse light-curves 
l|Marsh et al.ll2006l ) . Assuming the accretor obeys a fully de- 
generate WD M-R relation (a good approximation at the de- 
termined Mi), these authors determine Mi = 0.84±0.05A4q 




30 40 50 

P orb (minutes) 



Figure 19. The evolution of P or b for P or b > 25 min during the 
expansion phase. The lower panel shows P or b vs. P or b for two 
sets of Mtot each with two different log(i/> Cj i). Colours indicate 
log(-0 Cj i) and have the same meaning as in Fie. 1171 Dashed lines 
show systems with M\ t i = 0.35, M% , = 0.15Mq; short-dash dot- 
ted lines systems with M 1:i = 1.025, M%,i = 0.3M Q . The upper 
panel shows the relative number of systems expected in steady- 
state along each evolution track. As in Fig. 1171 the overall normal- 
ization is arbitrary, but the relative normalization between tracks 
is set by their P or b rates. The slowing of P or b during the donor's 
cooling phase leads to a peak in system numbers. The peak's P or b 
is diagnostic of Mtot (as well as n, although this is not show here). 
The magnitude of the peak increases with decreasing ip c ,i- 



and q = 0.035 ± 0.002, giving M 2 = 0.029 ± O.OO2M , 
~ 50% more massive than a zero-temperature WD that fills 
its Roche lobe at this system's P or b = 28.3 minutes. This M 2 
measurement provides us the first direct means of determin- 
ing an AM CVn donor's current entropy. From the deter- 
mined M2-range, we find that log(0 C) i) lies approximately 
in the range 1.60-1.35. If the actual range of He-star channel 
donor M-R evo l ution does not differ significantly from the 
iNelemans et al.l (|200ll ) fit for these systems, then a He-star 
channel system should have M2 « O.O5M0 at P rb = 28.3 
minutes. The evolve-MS channel also appears to produce a 
value of M2 which is too high for sy stems at this P OI .b (see Ta- 
ble 1 of lPodsiadlowski et alJ|20o3i . Thus SDSS J 0926+3624 
presents evidence that the WD channel indeed contributes to 
the observed AM CVn population and that this channel pro- 
duces non-zero entropy donors as predicted bv \Delove et all 

We can also predict the current M in SDSS 
J0926+3624. From our calculations, a system with the de- 
termined value of Mtot = 0.869M Q and M 2 = 0.029M Q at 
Porb = 28.3 minutes, has an M = 9.8 x 10~ 11 M Q yr" 1 . 
At fixed Porb when M 2 « Mi, M k M* /3 M% oc M? /3 q 2 
ilDelove et al.ll2005l ). so the error bars quoted in lMarsh et al.1 
fcoOrj ) for Mi and q provide an error of ~ 20% on M. Thus, 
we estimate M ~ 9.8 x 10" 11 ±2.0 x Kr u M yr" 1 . This M 
range is close to value at which the accretor's thermal evolu- 
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tion decouples from the compressio nal heating provided by 
the acc retion dBildsten et al.l l2006h . Thus, we can use Fig. 
1 from IBildsten et al.l 1 2006h to estimate the accretor's T c 
from M and find T c « 1.8-2.1 x 10 7 K. Using these M and 
T c ranges, we sum the accretor's cooling a nd compressional 
luminosity (see §2 of IBildsten et al1l2006l ) to estimate the 
accretor's T cB w 21,300-23,800 K, taking Mi = 0.84M Q to 
determine Ri ~ 0.01i?Q. More refined estimates of the ac- 
cretor's thermal properties will require more detailed calcu- 
lations taking into account variations in M evolution histo- 
ries, the time-dependent evolution of the accretor's envelope 
during the decoupling phase, and possibly the effect of He 
shell flashes on the accretor's surface. 

The measurements of Mtot and ip c ,i in additional sys- 
tems are required before the observational distributions of 
these parameters can be determined. However, the likelihood 
of discovering a system such as SDSS J0926+3624 given our 
theoretical models of the WD channel AM CVn population 
can be consider ed. We first compare M to t to the distribu- 
tion in Fig. 1 of lNelemans et all (|200ll ). The locus of points 
defining the region M to t = 0.869 ±O.O5M0 in this figure lies 
outside the M\-Mi parameter space these authors consider 
most likely. This raises the question of whether observing 
such a high M to t in SDSS J0926-3624 simply is the result 
of small number statistics or whether this is a hint of ad- 
ditional physics that skews the M to t distribution to higher 
values. An example of such physics would be the preferen- 
tial survival of high Mi,, s ystems during a di rect impact ac- 
cretion phase at contact ( Marsh et alj 12004! ) . Finally, from 



Figure El about 12% of WD channel AM CVn systems in 
our modelling have ip Ct i in the range consistent with SDSS 
J0926+3624, not taking into account the system's high Mtot • 
Given the rather flat f/> Cj i distribution expected from theory, 
determining M to t and ip c constraints in additional AM CVn 
systems is essential if we are to determine if this is consistent 
with our theoretical expectations or not. 

5.4 Observational Signatures of AM CVn Donors 

Here we consider predictions for the donor's contribution to 
the system's light in SDSS J0926+3624 based on our current 
models. The near edge-on inclination of SDSS J0926+3624 
makes this system a good candidate for discriminating the 
accretor's and donor's light without significant contamina- 
tion from the accr etion disk. The disco very spectrum of 
SDSS J0926+3624 (| Anderson et al.ll2005l ) is reproduced in 
Figure [20] from the SDSS archival data. It was obtained 
over a 3600s exposure — roughly two system orbits — and thus 
provides a phase averaged spectrum. It is decidedly non- 
blackbody in shape, indicating that a DB WD accretor is 
not the only contributor to the system's light. 

The question is whether the second component con- 
tributing to the flux is the donor or the disk? The 
iMarsh et all (|2006T ) Mi, M 2 constraints imply i?i w 0.017? Q 
and i?2 = O.O43i?0. From our current modelling, we pre- 
dict the donor's T e g ~ 1750-4500 K, correspon ding to an r\ 
range of 0.0-0.5. With i = 83.1° and q = 0.035 (|Marsh et all 
120061 ). the projected area of the disk's face ~ 2.77?2; i.e., out 
of eclipse roughly equivalent areas of disk and donor surfaces 
are seen. Thus, the disk may provi de a significant cont ribu- 
tion to the system's light. Indeed, IMarsh et al.l (|2006h find 
that the disk contributes ~ 50% of the flux in the r' band 
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Figure 20. Comparison between the Sloan spectrum for SDSS 
J0926+3624 (black line) and theoretical models. The upper panel 
shows a model spectrum (yellow line) calculated using contri- 
butions from a X e fj = 21, 000 K, R\ = O.OIR0 accretor (red 
line) and T off = 4400 K, R 2 = 0.043-Rq donor (which lies below 
the x-axis) at a distance of 293 pc. The lower panel shows the 
combined spectrum (yellow line) from a T e g = 39, 000 K accre- 
tor (red line) and a steady-state o-disk model (blue line) with 
M = 9.8 X -11 M Q yr _1 at a distance of 695 pc. 



and « 25% in the g' . They also find significant inter-orbit 
variability, so whether this particular ratio of flux contribu- 
tions is representative of conditions during the Sloan obser- 
vation is not certain. 

Given the prominent double-peaked He emission lines 
in this system, AM CVn phenomenology would argue the 
disk is either in a stable, seemingly optically thin, low-state 
or is an outbursting disk caught in quiescence. In either 
case, a stable, optically thick cv-disk model spectra is not 
expected. To check this, we calculated an a-disk spectra as- 
suming M = 9.8 x _11 Mq yr _1 and an outer radius of 0.7a. 
We added this to the accretor's flux modelled as a single 
T e ff blackbody. We then adjusted this T e £f and the system's 
distance to give a "by eye" best fit to the SDSS J0926+3624 
spectrum. The results are shown in the lower panel of Fig- 
ure !20l A rather hot, T c ft = 39, 000 K, accretor and a system 
distance of 695 pc is required for a reasonable fit. This model 
underestimates the continuum flux at wavelengths A < 4000 
Aand A > 8000 A. Additionally, the r' flux from the disk is 
significantly greate r than the accretor, in disagreement with 
IMarsh et all (|2006t ). This however may not be a significant 
issue given the source's variability. 

We also considered a model with accretor and donor 
contributions, but no disk. Both components were modelled 
as single T e ff blackbodies. The observed spectra is well fit 
by the combination of a T c ft = 21,000 K accretor and a 
T c ff = 4, 400 K donor at a distance of 293 pc (upper panel 
of Figure [20)) . T he accretor's T e g is a little below the range 
predicted by the IBildsten et al.l ([20061 ) theory (35T3J), while 
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the donor's T e g is consistent with an r\ somewhat below 
0.5. The fit in this case is somewhat better than the accre- 
tor+disk model. We should note that other acceptable fits 
can be found by collectively increasing or decreasing both 
components' temperatures and the system's distance. Thus, 
higher accretor T c g require higher r\ values for the donor to 
be hot enough to produce a good fit to the spectra. 

A parallax measurement for SDSS J0926+3624 will 
clearly distinguish between the stable a-disk+accretor and 
the accretor+donor model (as well as further constrain the 
component temperatures in this latter model). The better 
agreement between data and donor+accretor model and our 
expectation that the disk is not in a stable high-state already 
argues for the donor+accretor model. Given our poor under- 
standing of a low-state disk spectra and flux, however, one 
can only conclude that an accretor+donor model is fully 
consistent with both the observed spectrum and all theo- 
retical expectations for SDSS J0926+3624. Since this model 
doesn't include any disk contributions , it does fail to ex - 
plain the inferred disk properties of the lMarsh et al.l (|2006l ) 
results. It may be possible that it is a cool disk providing the 
long-A flux in this system, possibly exp laining both the Sloan 
spectrum and the lMarsh et al.l |2006l ) results. Better under- 
standing of quiescent He disks and phase-resolved spectral 
studies of this system may both be required to break this 
degeneracy between either a cold disk or donor as the source 
of long wavelength flux. 



6 SUMMARY 

We have implemented a new stellar evolution code in C++ 
that allows significant flexibility in use of input physics and 
defining systems of ODEs. This code has been used to model, 
for the first time, the full stellar structure of the donors 
in WD-channel AM CVn binaries. Specifically, the thermal 
and structural evolution of the donor has been calculated 
to determine how the donor's thermal evolution affects the 
evolution of the binary and to provide the first predictions 
for the donor's light in these systems. 

Systems forming through the WD channel are expected 
to have a range of donor properties, most i mportantly a 
range of initial (i.e., at contact) degeneracy (|Delove et al.l 
2005). We modelled the pr e-contact evoluti o n of W D chan- 
nel systems based on the iNelemans et al.l (|200lT ) popula- 
tion synthesis to determine suitable initial donor models 
for our subsequent calculations. The donors in this popu- 
lation have initial central degeneracy parameters between 
ip c<i Ril0-10 4 , with the distribution in f/) Ci , being rather flat 
between ~25-4000. Most of these systems make con- 
tact at P or b «2-ll minutes. This range of donor parameters 
is likely dependent on assumptions about the binary evolu- 
tion leading to their formation (in particular the CE-event 
prior to the AM CVn phase) and our assumptions about the 
donor's cooling during the pre-contact phase. Quantifying 
the potential effects of these assumptions will require more 
extensive modelling of the proto-donors in this pre-contact 
phase. 

Our evolutionary calculations show that WD-channel 
AM CVn systems have three phases of evolution. An M 
turn-on phase, during which M grows to its maximum value 
while R2 contracts under mass loss, and P or b decreases. This 



behaviour produces a turn-on phase lasting significantly 
longer than previous estimates (e.g., iMarsh fc Nelemans! 
l2005l ; IWillems fc Kalogerall2005h : ~ 10 4 -10 6 yrs depending 
on the donor's initial tp c ,i. In the second phase, the donor ex- 
pands adiabatically under mass loss. The third phase begins 
once the mass loss rate and the donor's thermal time have 
have both decreased enough to allow, starting at P or b ~ 45 
minutes, the donor to cool and contract to a fully degenerate 
configuration. We discussed how the system's initial condi- 
tions influence the later evolution of Ra(M2) and M(P or b). 
We also revised the upper limit to the M — P or b relation 
for WD-channel systems (given our initial condition deter- 
minations) . Finally, we predicted the donor's intrinsic L and 
T ft; during the adiabatic phase, L «10 -6 -10 -4 L@, while 
T efl fslOOO-1800 K. 

The flux generated by the accretion flow in these sys- 
tems can easily dominate the donor's intrinsic thermal out- 
put. We self-consistently modelled the impact of the accre- 
tion light on the donor's and binary's evolution in the grey 
approximation. The irradiation reduces the donor's intrin- 
sic L, producing a delay in the onset of cooling and slow- 
ing the donor's contraction once cooling does begin. This 
shifts the downturn in the M(P or b) relation seen in the non- 
irradiated donors to longer P or b- Irradiation also elevates 
the donor's photospheric temperature and surface luminos- 
ity (up to 10,000 K and 1O _2 L0, respectively, during the 
adiabatic phase). The observational signatures of irradiation 
depend on the irradiating flux's spectra, the opacity sources 
in the donors atmosphere and the efficiency of day-to-night 
side energy flow, all considerations beyond this work's scope. 
Our predictions for the irradiated donor's light are thus ap- 
proximate and valid in the limit of complete redistribution 
of flux in a grey atmosphere with a Bond albedo of zero. 

In comparison to prior predictions using sim- 
pler donor models, we find that previous assump- 
tions of fully-convective donors and adiab atic evolution 
(|Delove fc Bildstenl 120031 ; iDelove et all 120051 ) are not valid 
over much of the AM CVn evolution. The donors only be- 
come fully convective at P or b ~ 40 minutes and Af2 ~ 0.01. 
Prior to this, their shallower entropy profile leads to a slower 
expansion rate relative to fully convective models. Due to 
this, we argue that prior to the onset of cooling, M mea- 
surements could distinguish between systems formed in the 
H e-star versus WD c hannels (as opposed to the conclusions 
of lDelove et al.l [20051 ). After the cooling phase develops, all 
initial donor entropy information, including any that distin- 
guishes formation channels, will be erased. 

A system's evolution about its P or b-minimum and the 
occurrence of donor cooling phase both leave diagnostic 
signatures on the AM CVn population's P or b-distribution. 
Most significantly, the evolution of P or b — > at the Porb- 
minimum and the slowing of P or b evolution that occurs dur- 
ing the cooling phase lead to peaks in the P or b-distribution 
whose location and size depend on initial parameters (and r\ 
in the case of the cooling peaks). These could provide obser- 
vational diagnostics of the distribution of these parameters 
in the galactic AM CVn population. 

Finally, we showed that recent measurements of M\ and 
M 2 in the ec l ipsing AM CVn system, SDSS J0926+3624 
(|Marsh et alj |2006l ) provide direct evidence that WD- 
channel systems contribute to the AM CVn population and 
that this channel produces non-zero entropy donors as pre- 
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dieted bv lDelove et~ai1 <|2005h . We compared predictions for 
this system's light based on our models, showing that a com- 
posite spectrum consisting of donor and accretor contribu- 
tions is fully consistent with this system's discovery spectra. 
Based on this, we predict this system lies at a distance of 
« 290 pc. 

Current investigations include applying this theory to 
interpret the improving observational constraints in many 
of the known AM CVn systems (e.g., Roelofs et al. 2007, 
submitted) . In addition, we are considering how the physics 
relevant to the AM CVn system formation and early contact- 
phase survival could be probed by future LISA observations 
of the P or b-distribution of sources at P or b < 15 minutes. 
Other indicated work includes proper non-grey, phase de- 
pendent modelling of the donor's irradiated atmosphere as 
well as progress on understanding the contribution of low- 
state He accretion disks to the system's emission. 
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APPENDIX A: THE R 2 RESPONSE TO MASS 
LOSS 

The donor's R 2 can be expressed as the integral of equation 
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M = dm 
10 P 

The derivative of equation (|A1|) with respect to M 2 is 
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where Xp = {dhxP/d\sxp) T , xt = (<9 In P/<9 In T) p and V' = 
(d In T/d In P) m gives the change in T at fixed m that occurs 
due to changes in M 2 . 

The quantity {d\nP/dM 2 ) m ~ dlnP/dm ~ P C /(PM 2 ), 
showing the donor's outer layers dominate the R 2 response 
to M 2 . Also, (<91nT/<91nP) m > and V < V ad generically 
in our donors, so the local response to mass loss at fixed m is 
an expansion in r. The surface term t hen is the only d river of 
R decrease under mass loss (see also iHiellmingl [19891 ). Note 
that both these behaviours are apparent in Fig. [4] 

How the R 2 evolution depends on the donor's structure 



is more easily seen by transforming to P-coordinates and 
considering only the donor's outer layers. The contribution, 
SR 2 , to R 2 between surface at P p hot and some pressure p, ^> 
Pphot: 
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where gt ~ const, is the gravitational acceleration at Pf,. 
The change in this layer's thickness under mass loss is then: 
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(A4) 



where we have neglected for simplicity the surface 
term's contribution. The layer's response the depends on 
(din p/dM 2 )p, which can be rewritten as 
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APPENDIX B: THE LUMINOSITY PROFILE'S 
RESPONSE TO MASS LOSS 

To examine how the donor's /-profile evolution in response 
to mass loss depends on the donor's structure and relative 
ordering of r t h and r m , taking e = we can rewrite equation 
© as 
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where the first approximation holds when mass loss effects 
dominate (d In P/dt) m and the elapsed time under consider- 
ation St T t h. In the outer layers, the region most heavily 
weighted by the din P/dM 2 term (see Appendix E)) . this last 
expression can be expressed approximately as 
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Where St <g rth <S T m , heat transport is able to ap- 
proximately maintain the original thermal profile, produc- 
ing V' ~ V. For radiative regions, V < V a( j, producing 
d\nl/dm < as mass elements must absorb flux in order 
to increase their entropy as they move outward. In convec- 
tive regions, din I /dm > 0. Since the fractional change in 
I, SI 1 1 due to advection goes as ~ r t h/T" m , the overall flux 
decrement is rather small in this limit. 

In the opposite limit, when r t h 2> T m , mass elements 
are advected outward nearly adiabatically. Thus V' — > V a d 
from below in radiative regions and from above in convec- 
tive zones. Since r t h/T m can become very large, significant I 
perturbations due to mass loss can be driven by only a very 
slight non-adiabaticity in the advective flow. Our numerical 
calculations bear this out, showing that the magnitude of 
I perturbations are largest during the adiabatic mass loss 
phases. 
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